Dynamical Quantum Hysteresis 



OS 

o 
o 

(N 



(N 



Arnab Das 

Condense Matter and Statistical Physics Division, 
The Abdus Salam International Center for Theoretical Physics, Strada Costiera 11, 34014 Trieste, Italy 

(Dated: July 24, 2009) 

Here we report for the first time, the phenomenon of dynamical quantum hysteresis (DQH) and 
some interesting effects associated to it. We demonstrate the phenomenon in one-dimensional Ising 
model driven by a sinusoidal transverse field (symmetric about zero) at high frequency u> at zero 
temperature. DQH is manifested by the symmetry-breaking characterized by a non-zero value of 
the long-time average of the transverse magnetization, (the dynamical order parameter) close to its 
initial value. Interesting effects, namely, non-monotonic behavior of the order parameter with lo 
(showing curious peak structures), and stronger localization with increasing driving amplitude (in 
contrast to classical dynamical hysteresis) are observed. All the features are deduced analytically 
with exact asymptotic results. The analytical results are in excellent agreement with exact numerical 
results calculated for a system of 10 4 spins. Different experimental set-ups are discussed where the 
phenomenon may be observed. 



I 

C 

O 
o 



> 

(N 
O 

o 



X 



Hysteresis is a characteristic non-equilibrium response 
of a system to a periodic drive with the driving period 
much shorter than the intrinsic relaxation time of the sys- 
tem. Classically, a macroscopic relaxation time for the 
system to follow the instantaneous equilibrium state in a 
driven system often derives its origin from the presence 
of available metastable states (local energy /free-energy 
minima) that hold the system back as the field moves 
on. Hysteresis is hence observed, for example, when 
a system is driven about a first order phase transition 
point, [l], H, H, H|, or in presence of quenched disorders 
[5(. As a generic representative of the classical scenario 
let us consider a classical Ising model below the critical 
temperature Tc, in a longitudinal driving field oscillating 
symmetrically about zero. The drive will periodically el- 
evate and lower the two free-energy minima (correspond- 
ing to upward and downward magnetization), one with 
respect to the other. The quasi-static (hysteresis-free) 
dynamics would require the system following the global 
instantaneous free energy minimum through out the evo- 
lution. But since there is a free-energy barrier between 
the two minima, there is a macroscopic relaxation time 
for the system to relax from the higher minimum to the 
lower one (favored by the instantaneous hi), and if the 
period of oscillation is lower than this relaxation time, 
the system fails to follow the instantaneous field. For 
moderately high frequencies, there will be a phase-lag 
between the field and the instantaneous magnetization, 
giving rise to a symmetric hysteresis loop. But for very 
high frequencies, the system will be trapped dynamically 
in one of the minima, resulting in a non-zero value of the 
long-time average of the magnetization - a phenomenon 
known as dynamical (classical) hysteresis (DCH) 

Here we show for the first time that a similar dynami- 
cal localization may be observed in a periodically driven 
quantum system at zero temperature, due to radically 
different mechanisms (we call it dynamical quantum hys- 
teresis or DQH). Instead of an intrinsic relaxation time 



set by some available local energy-minima or metastable 
states, here the inherent timescale that competes with 
the driving period is set by the slowness of the evolu- 
tion required for the instantaneous state of the system to 
follow the ground state of the time-dependent Hamilto- 
nian adiabatically. Strong localization may occurs in the 
quantum case as a result of loss of adiabaticity, and addi- 
tionally, due to dynamics-dependent destructive interfer- 
ence effect, even in absence of any local energy minima 
and energy barriers. The phenomenon is completely dif- 
ferent from the quantum effects on hysteresis reported so 
far where in contrast to our case, quantum resonant 
tunneling actually helps relaxation. 

To demonstrate the phenomenon we consider one di- 
mensional transverse field Ising model (TFIM) [3] , which 
played a bench-marking role in hosting a large number of 
important studies in non-equilibrium quantum dynamics 
8]-[19], mostly concentrating around the effect of linear 
quench across the critical points/line on the density of 
defect produced, giving rise to the quantum versions of 
the Kibble-Zurek effect HUSH. The Hamiltonian for the 
TFIM with a sinusoidal transverse field may be given as 
follows. 
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(1) 



h z (t) = ho cos (uit). 



Here, a x J z represents the x/z component of the Pauli 
spin sitting at the z-th site, and h z (t) denotes the trans- 
verse field oscillating in time t with frequency u) and 
amplitude Hq. Hamiltonian |T]), with periodic boundary 
condition, when projected to the even-fcrmionic subspace 
(where lies the global ground state) via Jordan- Wigner 
transformation followed by Fourier transform, reduces to 
a direct sum of TV/2-number of 2 x 2 Hamiltonians H^'s 
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FIG. 1: Top panel (color online) (a) shows the variation of trans- 
verse magnetization m z (t) with time, for the driving amplitude 
ho = 10 (color online). The frequencies in decreasing order of av- 
erage m z (from top to bottom) are u> = 8.31661..., ui = 3.62305..., 
u) = 2.31114..., u) = 10.0, u> = 3 and u) = 5. Among them, 
the first three correspond to the first three zero's of Jo(2ho/to) 
at 2ho/w = 2.4048..., 5.5200... and 8.6537... respectively. The 
middle panel (b) shows (m z ) which is the long-time steady av- 
erage of m z , as a function of u) for ho = 10. Here the peaks 
Pi, P2 and P3 corresponds to the frequencies ui = 2.31114..., 
u) = 3.62305... and uj = 8.31661... respectively. Analytical results 
(m z ) = j (2fe /")l * S com P are d with exact numerical result for 
10 4 spins. The lower panel (c) shows the same as in (b), for ho = 20. 
The peaks Pi,P2 and P3 corresponds to u> = 3.392.., 4.622.., and 
7.246.. respectively, which corresponds to the zeros of Jo(2ho/u>) 
at 2h /u) = 11.971.., 8.653.., and 5.520... respectively. 



(see, for example, [HI]): 
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0, 1, ...,N/2 and^ fe = h z (t)+cosk 



and Afe = sinfc. Here Cj_ fc and c± k are respectively the 
creation and the annihilation operators for the Jordan- 
Wigner fermions (fc-fermions). In this representation, the 
time-dependent wave function IV'(i)) for the system may 
be expressed as a direct-product of the two-dimensional 
time-dependent wave functions: \ip{t)) = ®fc>o IV'fcW) 
in k space, with 

\Mt)) = Mfc(i)|0 +fe 0_ fc ) + v k (t)\ + k, -k) (3) 

where |0+fc0_fc) and | + k,—k) represents respectively 
the unoccupied and the doubly occupied states of ±/c- 



fermion. The time-dependent transverse magnetization 
in this notation reads 

m z {t) = W)\°*W)) = ^ E \ Vk ~ 1 ( 4 ) 

fc>o 

(dropping the site-index for translational invariance). 
Now we employ the unitary rotation \ip k (t)) = 
exp [|(2t cos k + sin Lot)a z ] \ip' k {t)) an d perform a 
simple expansion in terms of Bessel's function: 
exp [izsin6>] = Y,t=-oc J n{z)e m9 , J n (z) being the 
Bessel's function of the first kind with integer order n 
1221] . The rotated wave function 



\^' k (t))=<(t)\0+k,0-k)+v' k (t)\ + k,-k) 
follows the time-dependent Schrodinger equation 
. d 



l Jt Wk[t)) 



with the transformed Hamiltonian H' k (t) given by 
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(5) 



(6) 



(7) 



(8) 



We note that \u' k (t)\ 2 = \u k {t)\ 2 and \v' k {t)\ 2 = \v k (t)\ 2 . 

In order to investigate the localization, we start with 
the ground state of the Hamiltonian H(t = 0) in Eq. 
(fTJ) with ho 3> 1 (highly polarized in -f-z-direction). The 
dynamical localization may now be measured simply by 
monitoring the long-time average of m z : 



= lim — 

Tf^oo Tf 



m z (t)dt 



(9) 



The resulting time evolution of m z , obtained by exact nu- 
merical integration (for 10 4 spins) of the time-dependent 
Schrodinger equation for many sweeps as shown in Fig. 
[TJ panel (a) . The figure shows that that in each case, m z 
remains confined within a narrow range in the positive 
sector (starting with m z ~ 1) for all time, though the 
field h z oscillates symmetrically about zero. Moreover, 
its long-time average (m z ) is not a monotonic function 
of u>, rather, it shows some peaks appearing at certain 
values of the driving parameters (uj and h ) as shown in 
(Fig. □ b,c). 

This may be explained in terms of the dynamics of the 
fc— modes. Eq. of the form ([6]) can be solved analytically 
under resonance approximation (RA), which amounts to 
ignoring all the terms in the off-diagonal sum in H' k (t) ex- 
cept for the resonant term n — n r , for which the effective 
frequency £l k = \nui + 2cos/c| is the smallest. All other 
faster oscillating terms are neglected. With the initial 
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condition mentioned above (|«fc(0)| 2 ss 1), the solution of 
Eq. ([6]) for any n — n r under the RA gives [2^ : 



M*)l 2 = i- 



4Al 



AA% + (2cosfc + n r u>)'- 



■sin 2 Kt), (10) 



where = J„ r (2/i /w)Afc and vy. = 

A:A\ + (2cosfc + n r uj) 2 . In the concerned high 
frequency limit (ui ^> 2|cosfc|), we have n r — 0. From 
Eq. I|10p one sees (putting n r — 0) that the critical 
modes (k ~ 0, 7r), for which Aj, = sin(fc) ~ 0, |^fe(t)| 2 
oscillates with a vanishing amplitude proportional to 
A 2 and thus contribute to the freezing of m z around its 
initial value. However, for high u, the contribution from 
the off-critical modes (A& 3> 0) is nontrivial, and results 
in interesting behaviors as follows. For these modes, 
|wfe(t)| 2 oscillates with an amplitude A\j{A\ + cos 2 (A:)) 
about an average value 1 — A\j2{A\ + cos 2 (fc)). For 
Jo(2ho/uj) ^> 1, these modes follow the field almost 
adiabatically, but becomes more non-adiabatic as the 
the factor Jo(2ho/uj) becomes smaller. Taking the 
continuum limit of Eq. ([4J and integrating over k, one 
obtains a simple formula for (m 2 ): 





\ 1 / 
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l + |J (2/i /w)|. 



1(11) 
(12) 



This matches surprisingly well with the peaked- 
structured profile of (m z ) obtained by exact numerical 
integration, as shown in Fig. [T] (b) and (c). Here the 
time-average in the Tf — > oo limit is obtained analyti- 
cally by simply adding the time-average over a complete 
cycle for each mode, since Vk is finite for all k. The peaks 
occur for certain combinations of u> and ho, for which 
Jo (^ a ) = 0. Under this condition, the entire spectrum 
freezes, resulting in an absolute localization of the system 
around its initial state for all time. The condition for the 
dynamical freezing of a two-level system is well-known 
in the literature as the coherent destruction of tunnel- 



ing or CDT [24j . Such dynamics-dependent localizations 
may be interpreted as destructive interference of paths 
due to periodic driving, as discussed in [25| . In the limit 
lu — > oo, {m z ) — > | (due to the residual Rabi flopping 
dynamics). 

There is another interesting feature of DQH, which is 
also somewhat in contrary with the classical picture. In 
DCH, as one increases the strength of the driving field, 
it tends to enforce the system to follow it more closely, 
and for a strong enough driving field dynamical hysteresis 
cease to exit [4| . For example, a very strong longitudinal 
driving field in a classical Ising magnet would effectively 
flatten the free-energy barrier(s) between the metastable 
states created as a result of cooperative interaction, and 
would force the system to follow the field even without 
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FIG. 2: The plot shows variation of (m z ) (numerical) with hg 
for different cj. The counter-classical trend of higher degree of 
dynamical localization at higher driving amplitude is clear in the 
monotonically increasing nature of (m z ) with hg. 



much phase lag [26|. In the DCH studies, the localized 
phase is always observed below a certain value ho(w,T) 
of the amplitude of the driving field. As the amplitude 
is increased beyond that, the dynamical hysteresis is de- 
stroyed and the symmetric phase appears. 0. But in 
DQH, just the reverse trend is observed, as shown in Fig. 
[5J The figure shows that a higher freezing (characterized 
by higher values of (m z )) is obtained as ho is increased. 
In the high frequency limit, one can make the approx- 
imation Jq(x) w ^/^cos (a; — 7r/4) for x | and the 
expression (fT2|) reduces to 



(m z ) w (m z ) c 



where (m z ) oc = lirmj _> c 
localization in this limit. 



/UJ cos 



'2h 



'2h a 



(13) 



Eq. 



1, giving an absolute 
shows (m z )^r, — (m z ) 



vanishes as h 2 with a sinusoidal modulation as ho — > 
oo for a given uj. The classically intuitive picture of a 
stronger drive compelling the system more strongly to 
follow it, does not hold any more. 

In the high frequency limit considered, the contribu- 
tion of the off-critical mode to a non-zero (m z ) is high, 
since a vast part of the spectrum (even away from the 
critical modes) looses adiabaticity and freezes. As one 
goes to lower frequencies, the picture of freezing still 
continues to hold, since the non-adiabaticity of some de- 
gree is always guaranteed by the closing gap, but in this 
regime, DQH will be dominated mainly by the critical 
behavior of the system and peaks will be subdued (to be 
published later). 

Experimental verification of the DQH phenomenon can 
be achieved in several ways. First, realization of tunable 
transverse Ising model using trapped ions [27[ opens up 
a way of verifying our results directly. Similar realiza- 
tion is possible within the broad framework for realizing 
lattice-spin models with polar molecules in optical lat- 
tice (28| . In these systems, the exchange interaction J 
have experimental upper-limits (respectively of the order 
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of 22.1 kHz and between 10 — 100 kHz for the cases men- 
tioned above) but can be made arbitrarily small. Hence 
the range of high u> referred here (expressed in the units 
of J) may be brought down to any comfortable range, 
say of the order of few kHz in both the realizations. The 
phenomenon of DQH described here can be translated in 
two and three-dimensional Kitaev model 29|, 3(| , the ex- 
perimental realization of the former one being proposed 
already [3l| . The study of such hysteresis in random sys- 
tems in the context of quantum annealing 32|, l33l . l34j | 
with cyclic schedules would be interesting. The effect of 
DQH in quantum engineering, particularly, in ac-driven 
quantum systems [35| might be enormous. 

To summarize, we report the phenomenon of dynami- 
cal hysteresis in quantum system for the first time. The 
phenomenon exhibit some interesting features, namely, 
non-monotonic frequency dependence with interesting 
peak-structures, and completely reverse trend in ampli- 
tude dependence as observed in the classical cases. We 
demonstrate the basic features studying one-dimensional 
Ising model with oscillating transverse field. The demon- 
stration was entirely analytic and is in excellent agree- 
ment with the exact numerical calculations. 
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